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Abstract. The recent progress in the localized covariant density functional constrained by the 
relativistic Hartree-Fock theory is briefly presented by taking the Gamow-Teller resonance in 90 Zr 
as an example. It is shown that the constraints introduced by the Fock terms into the particle-hole 
residual interactions are straight forward and robust. 
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Introduction During the past decades, the covariant density functional theory 
(CDFT) has received wide attention due to its successful descriptions of both ground- 
state and excited-state properties of nuclei all over the nuclear chart. In this report, we 
will mainly focus on our recent progress in the localized covariant density functional 
constrained by the relativistic Hartree-Fock (RHF) theory [1]. 

Research on quantum mechanical many-body problems is essential in many fields 
of modern physics. By reducing the many-body problem formulated in terms of N- 
body wave functions to the one-body level with local density distributions p(r), the 
density functional theory (DFT) of Kohn and Sham [2] has accomplished great success. 
In particular, this theory states that there always exists a local single-particle potential 
Vks(r) so that the exact ground-state density of the interacting system can be reproduced 
by non-interacting particles moving in such local potential. 

The DFT has been also widely used in nuclear physics since the 1970s [3]. In par- 
ticular, the Lorentz symmetry taken into account in the CDFT puts stringent restric- 
tions on the number of parameters in the corresponding functionals without reducing 
the quality of the agreement with experimental data. During the past decades, the CDFT 
has been successful in describing a large variety of nuclear phenomena, including the 
equation of state in nuclear matter, ground-state properties of finite nuclei, collective ro- 
tational and vibrational excitations, as well as fission landscapes and low-lying spectra 
of transitional nuclei involving quantum phase transitions. The readers are referred to 
Refs. [4, 5, 6, 7, 8] for recent reviews. 

Part of these successes could only be achieved because the underlying functionals 
depend only on local densities, i.e., the relativistic Hartree (RH) scheme. However, a 
few common problems have been found in the isovector channels of the widely used 
RH parametrizations. Since a few years, the RHF calculations are possible for spherical 



nuclei all over the nuclear chart [9]. In contrast to the simple Hartree calculations, they 
reproduce successfully the proton-neutron effective mass splitting [9] (at the Hartree 
level it can be reproduced only if a iso vector-scalar meson is included [10]), shell 
structure evolutions [11, 12, 13, 14, 15], and spin-isospin resonances [16, 17, 18]. 
However, since the RHF theory introduces non-local potentials Vaj>(r,r'), its theoretical 
framework is much more involved and the simplicity of the Kohn-Sham scheme is 
lost. This restricts the applicability of deformed RHF theory to a few very light nuclei 
[19] and it prevents the treatment of effects beyond mean-field. Therefore, it is highly 
desirable to find a co variant density functional based on only local potentials, yet keeping 
the merits of the exchange terms. 

Recently, we started from the important observation that one of the successful and 
widely used RHF parametrizations PK02 [12] includes only three relatively heavy 
mesons, o, CO, and p, but no pions. Since the masses of these mesons are relatively 
heavy, the zero-range reduction for the corresponding finite-range Yukawa propagators 
becomes a reasonable approximation. As a step further, with the help of Fierz transfor- 
mation [20, 21], the Fock terms can be expressed as local Hartree terms. In such way, 
a Hartree-Fock equivalent covariant density functional can be derived while the con- 
straints introduced by the Fock terms of the RHF scheme can be taken into account. 
This RHF equivalent parametrization thus obtained is called PK02-H [1]. For details, 
the formalism for the effective Lagrangian density as well as the zero-range reduction 
and Fierz transformation can be found in Ref. [1]. 

In the following, we will take the Gamow-Teller resonance (GTR) in 90 Zr as an 
example to illustrate the ideas of the localized RHF equivalent RPA and the main 
conclusions. 

Results and discussion In the left panel of Fig. 1, the strength distribution of GTR 
in 90 Zr calculated by the particle-hole (ph) residual interactions of PK02-H is shown 
as solid line. It is found that, based on the same unperturbed single-particle spectra 
calculated by the RHF theory (dotted line), the present results by the localized RHF 
equivalent RPA calculations are quite similar to those by the original RHF+RPA [16] 
calculations (dashed line). The slight difference is due to the zero-range approximation. 
The experimental peak energy [22] is well reproduced. 

In order to understand the physical mechanisms, we present in the right panel 
of Fig. 1 the diagonal matrix elements of the most important GT ph configuration 
( Vi <?9/2) 1 ( n ^8i /2) 1 ■ The total matrix elements (full) are decomposed into the con- 
tributions from different channels. 

In the original RHF+RPA calculations with PK02, the ph residual interactions are 
provided by the o, CO, and p mesons. Both the direct and exchange terms of the ph 
residual interactions must be taken into account. It has been demonstrated that the total 
ph strengths are essentially determined by the delicate balance between the o and CO 
mesons [16, 18]. Since these two mesons are isoscalar, they can contribute via the 
Fock terms only, and correspondingly the effects of a and co mesons are repulsive and 
attractive, respectively. 

By taking the zero-range approximation of the finite-range meson-exchange Yukawa 
propagators and using the Fierz transformation, the ph interactions provided by the o, 
CO, and p mesons can be mapped into ten possible channels, i.e., the scalar (S), vector 
(V), tensor (T), pseudoscalar (PS), and pseudovector (PV) channels with the isoscalar 
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FIGURE 1. (Color online) Left panel: Strength distribution of GTR in 90 Zr calculated with the ph resid- 
ual interactions of PK02-H [1] with Lorentzian smearing parameter 1 MeV. The unperturbed excitations 
and the original RHF+RPA [16] results with PK02 [12] are shown as the dotted and dashed lines for 
comparison. The experimental peak energy [22] is denoted by the arrow. Right panel: Diagonal matrix 
elements of ph residual interactions for the (v 1^9/2) 1 {n^gi/iY configuration. The total matrix elements 
are decomposed into the contributions from different channels. 



and isovector (t) nature. For the charge-exchange modes, only the isovector channels 
contribute, since just the direct terms of the ph residual interactions are taken into 
account when the Fierz transformation has been preformed. 

Specifically, for the charge-exchange spin-flip modes like GTR, the dominant residual 
interactions are represented by the matrix elements of the spin-isospin operator [or] ■ 
[or] between the large components of Dirac spinors. They are composed of [c^t] ■ [oy/T] 
and [75 f] ■ [75 7 f] terms with the coefficients 

2a H=l4 and -c§ v = 14-74 + t4, (1) 

respectively, where a' 7 , 7, 7s are Dirac matrices, g x and m x (x = a, CO, p) are the meson- 
nucleon coupling strengths and meson masses in the original RHF parametrization 
PK02 [12], while a H are the Hartree equivalents. Correspondingly, the net contribution 
is then approximately proportional to 

2 «H-*~i(£j-i)' (2) 

As shown in the right panel of Fig. 1, it is true that, within the present localized RHF 
equivalent RPA approach, the total ph strengths are mainly composed of the positive tT 
and negative tPV contributions, while the contributions from the tS, tV, and tPS channels 
are negligible. In other words, the total ph strengths are determined by nothing but the 
delicate balance between the a and CO mesons (see Eq. (2)), which is well established 
in nuclear covariant density functionals. Therefore, the constraints introduced by the 
Fock terms of the RHF scheme into the ph residual interactions are straight forward and 
robust. 



Summary and perspectives In summary, a new local RHF equivalent covariant 
density functional was proposed recently in Ref. [1], where the constraints introduced 
by the Fock terms of the RHF scheme have been taken into account. In this way, the 
advantages of existing RH functionals can be maintained, while the problems in the 
isovector channel can be solved. This opens new and interesting perspectives for the 
development of nuclear local covariant density functionals with proper isoscalar and 
isovector properties. To follow this direction, the investigation on the relativistic local 
density approximation for the Coulomb Fock terms is in progress [23]. 
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